On the estimate of magnetic non-potentiality of sunspots derived 
using Hinode SOT/SP observations: Effect of polarimetric noise 

Sanjay Gosain, Sanjiv Kumar Tiwari and P. Venkatakrishrian 

Udaipur Solar Observatory, Physical Research Laboratory, Udaipur, Rajasthan- 313 001, 

O ! 

(N ■ INDIA 



3 

in 



6 



> 

in 

o 
in 

o 
o 



% 



Received ; accepted 



-2- 



ABSTRACT 



The accuracy of Milne-Eddington (ME) inversions, used to retrieve the mag- 
netic field vector, depends upon the signal-to-noise ratio (SNR) of the spectro- 
polarimetric observations. The SNR in real observations varies from pixel to 
pixel, therefore the accuracy of the field vector also varies over the map. The 
aim of this work is to study the effect of polarimetric noise on the inference of 
magnetic field vector and the magnetic non-potentiality of a real sunspot. To 
this end, we use Hinode SOT/SP vector magnetogram of a real sunspot NOAA 
10933 as an input to generate synthetic Stokes profiles under ME model assump- 
tions. We then add normally-distributed polarimetric noise of the level 0.5% of 
continuum intensity to these synthetic profiles and invert them again using ME 
code. This process is repeated 100 times with different realizations of noise. It 
is found that within most of the sunspot area (> 90% area) the spread in the (i) 
field strength is less than 8 Gauss, (ii) field inclination is less than 1 degree, and 
(iii) field azimuth is less than 5 degrees. Further, we determine the uncertainty 
in the magnetic non-potentiality of a sunspot as determined by the force-free pa- 
rameter a g and Spatially Averaged Signed Shear Angle (SASSA). It is found that 
for the sunspot studied here these parameters are a g = —3.5 ± 0.37 (xlO^m -1 ) 
and SASSA =—1.68 ± 0.014 °. This suggests that the SASSA is a less dispersion 
non-potentiality parameter as compared to a g . Further, we examine the effect of 
increasing noise levels viz. 0.01, 0.1, 0.5 and 1% of continuum intensity and find 
that SASSA is less vulnerable to noise as compared to a g parameter. 

Subject headings: Sunspot, Polarimetric noise, Magnetic field 
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Introduction 



Accurate determination of the vector magnetic field is very important for the studies 
of magnetic non-potentiality in active regions. The evolution of the active region magnetic 
field towards an increasingly non-potential state leads to buildup of magnetic free energy, 
i.e., energy above magnetic potential energy. This free magnetic energy is believed to drive 
the eruptive phenomena like flares and Coronal Mass Ejections (CMEs). Prediction of such 
phenomena on the Sun is very important for space weather forecasting, and requires good 
knowledge of the non-potentiality of the magnetic field in solar active regions. 

Force-free p arameter a g has been stu d ied to infer the non-po t entiality of sunspots 



for a long time (IPevtsqv et al. 



Hagino fc Sakurai 



( JHagino fc Sakurai 



2004, 



2005 



2004 ; 



1994 



Nandv 



Tiwari et al. 



995 



Abramenko et al. 



1996; 



Bao &: Zhang 



1998 



2006). Using the second moment of minimization 



2009aj) the relation between a g and the vector field 



components is given by the following relation : 



a,. 



{-^ \ fix nil ' i 



dy 



(1) 



The global alpha (a g ) actuall y gives twice t he deg ree of twist per unit axial length and not 
the twist (see Appendix A of iTiwari et al.l (J2009aJ)). Thus, if the length of the magnetic 
field structure in a volume is given, global twist can be obtained from alpha for force-free 
fields. 



Another parameter, spatially avera ged signed sh e ar ang le (SASSA), henceforth 



denoted as \1/, was recently proposed by 



Tiwari et al 



( J2009bl ) as a measure of magnetic 



non-potentiality in sunspots. This parameter is the spatial average of the angle between 
observed magnetic field and potential field azimuth. It is derived from the following relation. 



SASSA 



or 



* 



atan ( R R 4- R R ) ) 



(2) 



where B xo , B yo and B xp , B yp are observed and potential transverse components of the 
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magnetic field, respectively. The angled-braces represent the spatial average taken over all 
pixels except those below noise (see section 2.3). This parameter thus gives average shear 
angle on the photosph eric boundary and is independent of the force-free nature as well as 



shape of the sunspots (jVenkatakrishnan fc Tiwari 



20091 ). The high-quality Hinode data has 



also allowed these authors to study the contribution of local alpha v alues of umbr a 



penumbral structures to global alpha value of the sunspots. Further, iGosain et al. 



and 



(12009J ) 



evaluated \1/ of active region NOAA 10930 during X-class flare of 13 December 2006 using 



vector magnetogr ams obtained few 



(SP) instrument ( ILites et al 



2007 



(SOT) onboard Hinode satellite (IKosugi et al. 



lours before and afte r the flare by Spectro-Polarimeter 



Ichimoto et al. 



2007 



2008 ) with Solar Opti cal Telescope 



Tsuneta et al. 



20081 ). and found 



that the \1/ decreased after the flare. However, high-cadence vector magnetograms are 
needed to follow the evolution of non-potentiality, characterized by \1/, during the flare 
i nterval. Magnetogr aphs based on tunable filters like Imaging Vector Magnetograph (1V M) 



( IMickey et al. 



19961 ) and Solar Vector Magnetograph (SVM) (IGosain et al. 



2004 . 



20061 ) 



from ground and recently flown Helioseism ic and Magnetic Imager (HMI) on board Solar 



Dynamics Observatory (SDO) from space (jScherrer &: SDO/HMI Team 



20021 ) can provide 



such high-cadence magnetograms of active regions. 



Tiwari et al. 



( ]2009a[ ) had earlier evaluated the effect of polarimetric no ise on the 



1982). 



magnetic twist parameter a g and magnetic energy using a synthetic bipole (iLowl 
However, for real sunspots the distribution of magnetic field vector is not the same as in the 
case of synthetic bipole. Real sunspots have umbra and penumbra besides fine structure, 
as seen b y high resolutio n ground based and space based Hinode observations. In a recent 



study by 



Su et al. 



(120091 ) using Hinode observations it was found that the current and 
a distribution in a sunspot is not smooth but has a fine mesh-like structure with mixed 
polarity patches of a in the umbra and radial spine like pattern with alternating polarity 



in penumbra of a sunspot. At the umbra penumbra boundary they found an incomplete 
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annular ring with cu rrent and a values of sign opposite to that of global a i.e., a g . Further, 



Tiwari et al. 



(l2009bl ) have shown that the distribution of a in penumbra almost cancels 
to zero while in umbra there is a net value which bears same sign as a g and has same 
magnitude a s a ? . F urther, the new parameter \l/ that has been introduced recently by 



Tiwari et al. 



(]2009bl ) needs to be assessed for robustness in the presence of uncertainties in 



th e magnetic fi e 



of 



Tiwari et al. 



d para meters of a real sunspot. This work therefore extends the results 



( J2009al ) to the case of a real sunspot observed by Hinode and to the new 



parameter \l/. 

The uncertainties in vector magnetic field obtained by fitting the observed Stokes 
profiles with model profiles can arise mainly due to the following reasons : (i) invalidity 
of the model atmosphere used, e.g., ME model cannot fit asymmetric Stokes profiles 
which arise due to variation of physical parameters within the line forming region, and 
(ii) polarimetric noise, e.g., sensitivity in determining the field parameters depends on the 
Sig nal to Noise Ratio (SNR ) of th e Stokes profiles. The former reason has been evaluated 



by 



Westendorp Plaza et al. 



on Response-Functions; 



(200 1) by comparing ME and SI R (Stokes Inversion based 



Ruiz Cobo fc del Toro Iniestal ( 119921 )) inversions. SIR inversions 



yield variation of physical parameters within the line forming region, by exploiting the 
asymmetric nature of Stokes profiles. They have shown that using ME inverted field 
parameters is equivalent to using field parameters averaged over the line forming region 
obtained by SIR inversions. So, neglecting effect of Stokes asymmetry we focused on the 
errors in the field measurements due to polarimetric noise. 

The paper is organized as follows. In section 2 we describe the Hinode SOT/SP 
observations and the method of simulating noisy Stokes profiles and their inversion. The 
results are presented in section 3. The importance of the non-potentiality parameters 
and the inaccuracy in their determination as a result of polarimetric noise is discussed in 
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section 4. Section 5 gives the conclusions based on the present study and the future studies 
required. 

2. Data and Method 

2.1. Observational Data 

The high resolution spectro-polarimetric observations were made by the Spectro- 
Polarimeter (SP) instrument working with Solar Optical Telescope (SOT) onboard Hinode 
space mission. We choose a regular isolated sunspot in NOAA 10933 during 05 January 
2007, when it was located close to the disk center (// = 0.99) during 07:20 UT. The 
sunspot was scanned with "fast-map" observing mode by the SP instrument. In this 
mode the spatial sampling along the slit is 0.32" per pixel and across the slit is 0.29" 
per pixel with integration time of 1.6 seconds. The spectral lines used for polarimetric 
measurement are Fe I 6301.5 and 6302.5 A line pair. The vector magnetic field for 



this sunspot was derived using Milne-Eddington code nam ed MERLI1N 



, prov ided under 



2003). It is an 



Community Spectro-polarimetry Analysis Center (CSAC) ( ILites et al. 
inversion code based on the least-squares fitting of the observed Stokes profiles using the 
Levenberg-Marquardt algorithm and is quick due to parallel computing. It assumes a 
standard Milne-Eddington atmosphere to retrieve the magnetic field vector, line-of-sight 
velocity, source function, Doppler broadening as well as the macro-turbulence and the stray 
light filling factor. For Monte-Carlo estimation of erro r in magnetic pa rameters we use 



20041 ) . We compared 



another Milne-Eddington inversion code named Helix (jLagg et al. 
the magnetic field parameters retrieved by the two codes after inverting the SOT/SP 
dataset used in the present study. We made difference maps of the retrieved magnetic 
field parameters i.e., 5M = M MERLIN — M HELIX , where M = B (field strength), 7 
(inclination) or -0 (azimuth). Within the sunspot, the two codes retrieve identical magnetic 
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field parameters giving zero mean of SM and (1-cr) standard-deviation of 50 Gauss, 0.75 
degrees and 2.5 degrees in field-strength, inclination and azimuth angle, respectively. These 
values are comparable to the standard errors of inversion by MERLIN, reported in figure 
3. Standard errors mean that if same inversion code is repeatedly applied to a given 
observed (noisy) profile, it will lead to same model atmosphere within the error bars given 
by the standard error est imate, provided the minimum of the merit function is reached 



(IBellot Rubio et al. 



20001 ). So, the comparison suggests that the two codes are identical such 
that the differences in the output are of same magnitude as the difference of a particular 
code with itself. The reason for choosing Helix over MERLIN in our Monte-Carlo error 
estimation is the in-built routines for adding noise to the synthetic Stokes profiles in Helix 
code which helps us save time as we avoid adding noise offline. The azimuthal ambiguity 
of the tr ansverse comp onent of the magnetic field vector was resolved using acute angle 



method (IHarvey 



19691 ) 



2.2. The Monte-Carlo Approach for Determination of Uncertainties 

The uncertainties in the field parameters derived using ME inversion of observed Stokes 

profiles, purely as a result of polarimetric noise present in the observed data, are considered 

here. There are two ways of determining the uncertainties: 

(i) Standard Errors : Basically, in the Stokes profile inversions, a merit function is defined 

as 

4 m 2 

fe=l 4=1 kt 

where k—1 to 4 represent the four Stokes parameters, i— 1 to M represent the number of 
wavelength samples for observed and synthetic Stokes profiles, v represents the number of 
degrees of freedom (i.e., the number of observables minus the number of free parameters), 
and the vector a consists the parameters of the model, er^ represent the noise in the 



observations and Wki represents the weight given to the data points. This merit function, x 2 , 
is then minimized by using non-linear least squares method to obtain Xmin- The standard 
errors in the model parameters basically depend upon the curvatu re of the x 2 function near 



the region of minimum in the parameter space (IPress et al. 



19861 ). More details about the 



estimati on of the standard e rrors in the model parameters derived by Stokes inversion is 



given in 



Bellot Rubio et al. 



(120001 ). 



(ii) Monte-Carlo Errors : In Monte-Carlo method we create several artificial realizations 
of the spectro-polarimetric data. This is done by generating the synthetic Stokes profiles 
and adding different realizations of noise to them. These artificial datasets are then inverted 
to obtain the magnetic parameters, leading to a collection of results. The errors are then 
determined by finding the spread or standard deviation in the derived parameters. It may 
be noted that the Monte-Carlo errors have no obvious relationship with the standard errors. 



Westendorp Plaza et al. 



(120011 ) has shown that in case of SIR inversion of a single 
profile, the standard errors are slightly larger than the estimates from Monte-Carlo method. 
However, this depends on the distribution of field vector, and so the comparison of errors 
by the two methods may vary from pixel to pixel. Further, the Monte-Carlo error estimates 
depend upon the noise in the data as well as the noise in the process of inversion. Thus 
these are useful to establish the statistical significance of the results. 

We use Monte-Carlo approach to determine the statistical spread or uncertainty in the 
magnetic field parameters as well as the non-potentiality parameters of active region arising 
due to polarimetric noise. The approach is as follows: 

(1) First we take real observations of a sunspot by Hinode SOT/SP and invert it with a 
Milne-Eddington code to retrieve magnetic field vector. 

(2) resolve the azimuthal ambiguity in the transverse field using acute angle method. 

(3) consider this vector field as the "true field" [B , 70, ipo], and the derived quantities 



-9- 

a(i?o,7o,^o) and ^(B ,j ,ip ) as "true non-potentiality parameters" [ao,^o], for reference. 
(4) using the "true field" value at each pixel we generate synthetic Stokes profiles and add 
polarimetric noise of level 0.5% of continuum intensity. These profiles are then inverted 
again under ME approximation to derive field vector [B, 7, ip] and [a g , \P]. This process 
is repeated 100 times, each time with different realization of polarimetric noise. The 
uncertainty in the field strength, inclination and azimuth for each pixel and the uncertainty 
in the parameters a g and \l/ for the vector map is then estimated from the 1-a standard 
deviation in the 100 values of these parameters. This standard deviation is called the 
Monte-Carlo error. 

This simulation gives us an idea about the spread (standard deviation) in derived field 
vector and the non-potentiality parameters of a real sunspot, arising due to polarimetric 
noise. This knowledge of standard deviation in the magnetic non-potentiality parameters 
a g and \1/ is important to determine whether the observed changes in these parameters, for 
example in relation to flares, are significant or not. Such studies are yet to be done and 
the present work will establish the level of uncertainties in the parameter \l/ or a g due to 
polarimetric noise in modern observations, such as from Hinode. 

2.3. Method of Adding the Noise 

We first generate synthetic Stokes profiles corresponding to the "true field " using 



Milne-Eddington based Stokes profile synthesis and inversion code named Helix (ILagg et al. 



20041 ). To these synthetic profiles we add normally distributed random noise with the 3-cr 
level of 0.5% of continuum intensity (J c ). The noise, N(X), is added to the synthetic Stokes 
profiles, S syn (X), as follows: 
(1) First a pseudo random-number sequence is generated which is normally distributed with 
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a zero mean and a 3-cr standard deviation of given level, say L. In our case L=0.5% of I c . 

N(X) = Random_Number (Seed, L) 
(2) This sequence is then added to the synthetic Stokes profiles to yield noisy Stokes profile, 

&noi\X). 

S noi (X) = S syn (X) + N(X) 

The noise level is estimated from the observed signals in the continuum of the Stokes 
spectra. A continuum window between 6302.83 A and 6303.28 Ais selected for monitoring 
the noise in Stokes signal. The left panel of the figure [T] shows the histogram of noise in the 
observed Stokes profiles for a large number of pixels (1024 pixels) for a typical SOT/SP 
scan in "fast-map" observing mode. A Gaussian fit to this distribution yields a 3-cr value of 
0.5% of I c . This is the noise level that we used in our simulations. The histogram in the 
right panel of figure Q] shows the distribution of artificial noise that we add to the synthetic 
profiles. These profiles are then inverted with the Helix code. This process is repeated 100 
times and so for each pixel we have 100 values distributed around a mean value. 

While modern spectro-polarimetric observations have typically noise levels of the order 
of 0.5% of J c , archived observations from the ground-based instruments, which might be 
used for synoptic studies, may have higher levels of noise. Therefore, we also carried out 
an exercise to check the variation in the magnetic as well as non-potentiality parameters, 
\& and otg, with increasing polarimetric noise levels. We added four different levels of noise, 
i.e., 0.01, 0.1, 0.5 and 1 % of continuum intensity. 

2.4. Method of computing a g and \P 

The computation of the non-potentiality parameter a g and \1/ is done as proposed by 
Tiwari et al. (2009a, 2009b). The expression for computation of these parameters is given 
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in equations (1) and (2) in the introduction section. Only those pixels which have field 
values above a certain noise level are analyzed. This filtering is done in the following way: 
we select a quiet region on the Sun and evaluate 1-er deviation in the three vector field 
components B x , B y and B z separately. The box selected for this estimation is shown in 
top-left panel of figure [2 For transverse vector fields, we take 1-a noise level as the resultant 
deviations obtained in B x and B y . Only those pixels where transverse and line-of-sight 
(LOS) field both are together greater than twice the above mentioned noise level of \-a are 
analyzed. The \-a value for longitudinal and transverse field inside the box is 380 and 258 
Gauss, respectively. 

3. Results 

3.1. Effect of noise on magnetic field parameters 

The figure [2] shows the effect of noise on magnetic field parameters for the NOAA 
10933 sunspot. The top row shows the initial magnetic field parameters [B , 7 , ^o] derived 
by ME inversion of Stokes profiles obtained by Hinode SOT/SP scan during 07:20 UT on 
Jan 5, 2009. The middle row shows the map of Monte-Carlo error, i.e., the 1-<t standard 
deviation in the magnetic field parameters, derived using the Monte-Carlo method. The 
bottom row shows the map of standard error in the field parameters from the least-squares 
fit of the Stokes profiles. The maximum error (taking into account both the standard error 
and Monte-Carlo error estimate) inside the sunspot is less than 50 Gauss for the field 
strength. While for inclination and azimuth it is less than few degrees. The azimuth errors 
are typically largest in the umbral and plage region where the field is almost vertical and 
azimuth is not well defined. 

In figure [2J we isolate the umbral and penumbral regions by using continuum intensity 
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thresholds. We compare the standard errors and the Monte-Carlo errors in these two 

regions in the figure [3J The number of pixels in the umbral and penumbral regions plotted 

in figure [3] are 1953 and 13615, respectively. It may be seen that : 

(i) in the umbral and penumbral regions the Monte- Carlo errors are typically less than 5 

Gauss and 8 Gauss for field strength, less than 1 and 0.5 degrees for field inclination, and 

less than 5 and 2 degrees for field azimuth, respectively. 

(ii) the standard error is larger than Monte-Carlo er ror estimates in both the umbr al and 



Westendorp Plaza et al 



(120011 ) in case 



penumbral regions for field strength as was found by 

of the inversion of a single profile. While for field azimuth and inclination it has an opposite 

relation, except for middle panel of bottom row. 

The standard error depends upon the sensitivity of the \ 2 to the changes in a particular 
parameter. It can so happen that the Stokes profiles might not have significant sensitivity 
to the changes in that parameter. This results in a large standard error. On the other 
hand, the Monte-Carlo method examines the changes in the derived parameters resulting 
from the changes in the Stokes profiles arising due to random fluctuations in the profile as a 
result of different realizations of noise being added. Since both errors depend upon different 
origins of the profile fluctuations, they need not produce the same results. 

In general, considering all the panels in figure 3, it may be concluded that the standard 
errors are larger than the Monte-Carlo errors. Here, one should bear in mind that standard 
errors shown in figure 3, are for inversion of real Stokes profiles which possess asymmetry, 
apart from the polarimetric noise. While, Monte-Carlo error correspond to repeated 
inversion of purely synthetic, and therefore symmetric (or antisymmetric for Stokes-V) 
profiles with polarimetric noise. Therefore, the main source of difference between the two 
type of errors could be attributed to the presence of Stokes asymmetry in the real stokes 
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profiles leading to large value of standard error. Ideally, one would like to perform a 
Monte-Carlo error estimate by obtaining several simultaneous observational datasets of a 
real sunspot where the Stokes profiles would have inherent Stokes asymmetry as well as 
different realization of polarimetric noise. However, in the absence of such a possibility 
we used present method to study the effect of polarimetric noise alone, neglecting Stokes 
asymmetry. The results of the present study are therefore only the lower limits of the errors 
that can occur practically in real observations as well as in the non-potentiality parameters 
^ and a,g. 

3.2. Effect of noise on Azimuthal Ambiguity Resolution 

It is well known that the Zeeman effect diagnostics cannot detect the direction of the 
transverse field component and so a 180 degree ambiguity remains in the determination 
of the field azimuth. Various methods, however, have been developed by the researchers 
to resolve this ambiguity using different arguments (Leka et al. 2009). One of the 
most common and widely applicable method is the so-called acute angle method. 
In this method the angle between the observed and potential transverse field, i.e., 
9 = acos(B° bs .Bf° )/(| B° bs \\ Bf° |), is computed and the solution for which the value 
of 9 is an acute angle is considered as the correct solution. However, in the presence of 
polarimetric noise in the observations, this method also can fail specially in situations where 
the observed field is highly sheared. Such highly sheared regions are found near polarity 
inversion line (PIL) of the active regions. 

Fortunately, in a normal round sunspot, like the one used in our simulations, there 
are no high-shear regions and therefore the azimuthal ambiguity is easily solved with acute 
angle method. Nevertheless, we need to check the effect of polarimetric noise on the azimuth 
ambiguity resolution in our simulations before we examine the uncertainty in the \l/ and a g 
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parameter of the sunspot in 100 realizations. After resolving the azimuthal ambiguity using 
acute angle method for the 100 realizations of the sunspot vector maps, we made a map of 
1-er standard deviation for azimuth angle, as shown in middle-right panel of figure [2J If the 
ambiguity is not resolved properly in these 100 realizations, there will be fluctuations of the 
order of 180 degrees, leading to a large value of standard deviation in azimuth value for a 
given pixel. However, in the right panel of the middle row in figure [2] we see that: 
(i) the resolution of the azimuth angle is quite stable for the most part within the sunspot, 
specially in the penumbral region. Outside the sunspot, in the quiet and facular areas, the 
errors are large. This is mainly due to poor SNR in these areas as a result of weaker and/or 
vertical fields. 

(ii) As we go towards the umbra the values of standard deviation are large. We examined 
the values of the field azimuth in these regions and found that values do not vary to the 
extent of 180 degrees and therefore the large errors in umbral region is not due to ambiguity 
solver but due to poor SNR in Stokes-Q and U observations. Only in the very central 
part of the umbra, where the field inclination is close to 90 degrees, the azimuth loses its 
meaning and so we see a large spread. 

Further, we also checked the effect of ambiguity solver with increasing noise i.e., four 
different levels of noise of 0.01, 0.1, 0.5 and 1 % of continuum intensity. Here also the 
azimuth ambiguity resolution is stable for most part of the sunspot. 



3.3. Effect of noise on ^ and a g 

The figure H] shows histogram of \l/ and a g corresponding to 100 Monte-Carlo 
realizations of vector field. It can be noticed that: 

1. The \l/ is not affected much by noise: The distribution of \& corresponding to 100 
realizations show less scatter (~ 1%) with \I/ =—1.68 ± 0.014 °. 
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2. The values of a g are more affected: The distribution shows large scatter (~ 10%) in 
values with a g = -3.5 ± 0.37 (xlO^m- 1 ). 

These results show that the \& may be more reliable as compared to a g for typical 
polarimetric noise present in modern spectro-polarimetric observations as it is a lower 
dispersion parameter. 

We also check the effect of increasing noise in the polarimetric measurements on the 
magnetic as well as non-potentiality parameters of the same sunspot. In order to compare 
the effect of increasing noise on the field parameters in the sunspot we isolate the umbral 
and penumbral region, as was done in figure [2J and do a scatter plot between the input 
and output field parameters, for different levels of noise. The scatter plot for umbral and 
penumbral region is shown in figure [5] and figure [6] respectively. It can be seen that, in the 
umbra, where the field is mostly vertical, the azimuth determination is more affected with 
increasing noise than the field strength and inclination. Specially, for 0.5% and 1% noise 
levels the spread is large in the azimuth and inclination values. In comparison, the scatter 
plot for the penumbral region, where the field is not so vertical, shows that the azimuth is 
determined with less spread even for 0.5% and 1% noise levels, respectively. 

Table [1] lists the non-potentiality parameters, derived from the ME inverted magnetic 
field vector, after adding noise of different levels in the synthetic Stokes profiles. In Table 1, 
the spatial fluctuations of the shear angle \& as well as local a in the map, are represented 
as cry and a a . Starting with no noise and then adding random noise of 0.01, 0.1, 0.5 and 
1 % of continuum intensity in the Stokes profiles we get the following results: (i) The sign 
of the twist in the sunspot magnetic field, as inferred by both a g and \1/, is negative and 
is not affected even when the SNR of Stokes profiles is poor, (ii) the absolute value of \1/ 
tends to decrease systematically with increasing noise in the Stokes profiles, and (iii) the 
fluctuation of 2.5% in the \1/ = —1.65 ± 0.04° is less than the fluctuation of about 5% in 
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9^,-^ 



a g = 3.15±0.17(xl0- 9 m 

4. Discussion 

Developing quantitative measures of magnetic non-potentiality in solar active regions 
is very important for flare research. The line-of-sight magnetic field alone is insufficient 
for this purpose, therefore people started measuring vector magnetic fields. The stresses 
in the magnetic field were quantified in terms of the magnetic shear angle, by comparing 
the observed field with the potential- field (or current-free fields). The distribution of shear 
angle over active region showed that the filament bearing PIL regions are characterized by 
high-shear and are potential sites for flares. Shear, as a local measure of non-potentiality 
can be easily extended as a whole active region measure, like SASSA ^ (Tiwari et al. 
2009b). 

Other physical measures of whole active region non-potentiality are the virial estimate 
of free-energy and a g . Both of these measures require the photospheric magnetic field 
to be force-free. However, force-free condition may not be justified as the plasma (3 in 
the photosphere is not much smaller than unity and so the non-magnetic forces are not 
negligible. The effect of the polarimetric noise on the virial free-energy estimate was 
determined by Klimchuck et al. (1992) and Tiwari et al. (2009 a). Howe ver, the results of 



19821). While, in real 



these studies are valid for analytic force-free field solution of ( iLowl 
sunspots, the effect of non-force-freeness of photospheric field on virial energy and a g may 
be larger than the effect of polarimetric errors quantified by Klimchuck et al. (1992) and 
Tiwari et al. (2009a). The whole active region non-potentiality characterized by SASSA \l/ 
is, however, free of force-free assumption and in the present work we studied the effect of 
polarimetric noise on this parameter. 
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We simulated polarimetric profiles with noise level of 0.5% of I c , which is typical of 
Hinode SOT/SP observations as shown in figure 1. It must be noted that the synthesis of 
these profiles is done under the assumption of ME model atmosphere and so the present 
results pertain to symmetric Stokes profiles only. The asymmetry of the Stokes profiles 
would lead to a systematic rather than random effect which needs to be quantified in a 
separate study. In hundred realizations of such noisy symmetric Stokes profiles we found 
that the parameter \1/ is statistically more stable than a g . The reason for this stability of 
^ as compared to a g can be understood as follows. In Eqn. (2) \I/ is derived as a simple 
summation of angles. So, the random errors would cancel each other in summation. While, 
in Eqn. (1) we note that a g depends upon three components B x , B y and B z and is not 
a simple summation and therefore the random errors in the three parameters would not 
vanish statistically. Thus, a g would be more vulnerable to the noise as compared to \1/. 

Further, in the presence of different noise levels i.e., 0.01, 0.1, 0.5 and 1% of J c , it is 
found that the standard deviation in \1/ is less vulnerable to increasing noise than a g . The 
absolute value of \1/ tends to reduce systematically with increasing noise. 

Further, one should keep in mind that the accuracy of derived magnetic field parameters 
also depends upon the number of free parameters (parameters of the model atmosphere) 
used in the fitting procedure. The present results pertain to a single-component model 
atmosphere with stray-light component. The results would be different when one employs 
more number of free parameters, like two-component model atmosphere (Leka 2001; Lites 
et al. 2002) or depth- dependent (stratified) model atmosphere used in SIR inversions. In 
general the uncertainty in the best-fit model parameters will deteriorate with increasing 
number of model (or free) parameters. 
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5. Conclusion 

We generated an ensemble of artificial spectro-polarimetric dataset for a given sunspot. 
This ensemble of data is then inverted to give an ensemble of vector maps. These 
maps are then used to estimate of uncertainties in the field parameters as well as in the 
non-potentiality parameters, a g and \I/, of the sunspot. The stan dard errors of ME fittin g 



are gi ven by the inversion codes according to 



Press et al. 



(119861 ) and 



Bellot Rubio et al. 



( 120001 ) . The Monte-Carlo errors give us an independent method to cross-check the standard 
error estimates. In figure 3, we show that the Monte-Carlo errors for field strength, 
inclination and azimuth in the sunspot is determined within ± 8 G, ± 1 degree and ± 
5 degree, respectively. A comparison shows that the Monte-Carlo errors in field strength 
are typically smaller than the standard errors, while it is opposite for the field azimuth 
and inclination, except for the field inclination in penumbra where again the Monte-Carlo 
errors are less than standard errors. In general, the standard errors are more conservative 
estimates because they include the effects of polarimetric noise as well as those of Stokes 
asymmetry, while Monte-Carlo errors account only for the polarimetric noise. 

The effect of polarimetric noise on the parameters characterizing the non-potentiality 
of a sunspot magnetic field suggests that \1/ is more robust than a g (as shown by the 
histogram in figure 4). Further, \1/ appears to be a stable parameter with increasing noise 
in the polarimetric data (Table 1). 

Thus, ty will be useful for studying the evolution of non-potentiality of the active 
region magnetic fields, which in turn, can help in the prediction of flare occurrence. The 
parameter \1/ in a large number of active regions as well as the evolution of \1/ in a flaring 
regions will be evaluated in a future work. 

We thank the referee for detailed comments and suggestions. We thank Dr. Bruce 
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as international partners. It is operated by these agencies in co-operation with ESA and 
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Fig. L— The left panel shows the histogram of the noise in the Stokes profiles (S/I c ) 
for typical SOT/SP "fast-mode" scan. The right panel shows the histogram of normally 
distributed psuedo-noise added in the synthetic profiles for Monte-Carlo simulations. The 
solid line shows the fitted Gaussian with 1-er noise level is about 0.17% of I c . 
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Fig. 2. — The top row shows the maps of magnetic field strength, inclination and azimuth for 
the sunspot in NOAA 10933 derived using Milne-Eddington inversion. The solid and dashed 
line contours mark the penumbral and umbral regions, respectively. The white box in top 
left panel marks the box used for estimation of noise levels in magnetic field components. 
The middle row shows the error maps (1-cr standard deviation) for field strength, inclination 
and azimuth derived using Monte-Carlo method. The bottom row shows the maps of the 
standard error in Milne-Eddington fitting for the field strength, inclination and azimuth. 
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Fig. 3. — Top panel: The scatter plot between the Monte-Carlo errors and standard errors 
of fitting are shown for field strength, inclination and azimuth for all pixels in umbral region. 
Bottom panel: The scatter plots for pixels inside penumbral region. 
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Fig. 4.— The histograms show distribution of \& and a 9 corresponding to NOAA 10933. 
The vector maps were derived hundred times, and each time different realization of normally 
distributed noise (3-cr = 0.5% of I c ) is added to the synthetic profiles. 
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Fig. 5. — The scatter plots between the input and derived magnetic field parameters after 
different level of noise is added to Stokes profiles. The top of each panel shows the 3-cr level 
of normally distributed noise that is added to the synthetic profiles. The x-axis shows the 
input magnetic field parameter used in the profile synthesis and the y-axis shows the output 
value of the parameter derived from the noise added synthetic profile. 
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Fig. 6. — Same as figure 7, but for penumbral region. 
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Table 1: Effect of increasing polarimetric noise in the estimation of \1/ and a 6 

Noise (% of I c ) $ (°) <j 9 (°) a^m" 1 ) ^(m" 1 ) 

-1.692 15.724 -3.334 x 1(T 9 3.572 x 1(T 7 

0.01 -1.688 15.726 -3.360 x 10~ 9 3.573 x 10~ 7 

-1.685 15.702 -3.221 x 10~ 9 3.574 x 10~ 7 

-1.648 15.636 -2.976 x 10~ 9 3.562 x 10~ 7 

-1.604 15.364 -3.026 x 10~ 9 3.539 x 10~ 7 
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